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Recent experiments have shown that it is possible to create an in-plane harmonic 
potential trap for a two-dimensional (2D) gas of exciton-polaritons in a microcavity 
structure, and evidence has been reported of Bose-Einstein condensation of polari- 
tons accumulated in this type of trap. We present here the theory of Bose-Einstein 
condensation (BEC) and superfluidity of the exciton polaritons in a harmonic poten- 
tial trap. Along the way, we determine a general method for defining the superfluid 
fraction in a 2D trap, in terms of angular momentum representation. We show that 
in the continuum limit, as the trap becomes shallower the superfluid fraction ap- 
proaches the 2D Kosterlitz-Thouless limit, while the condensate fraction approaches 
zero, as expected. 

PACS numbers: 71.36,+c, 03.75.Hh, 73.20.Mf, 73.21.Fg 



I. INTRODUCTION 

In the past decade, there has been extensive work on Bose coherent effects of two- 
dimensional (2D) exciton-polaritons in cavitie. (For general reviews see Refs. El and Q). 
A microcavity is formed by two mirrors opposite each other, as in a laser cavity, with 
quantum wells embedded in the cavity at the antinodes of the confined optical mode. The 
resonant exciton-photon coupling leads to two polariton branches in the spectrum. The 
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lower polariton (LP) branch has a minimum at k = with a very small effective mass, in 
the range 10~ 5 — f 0~ 4 of the vacuum electron mass, depending on the details of the structure. 
These quasiparticles act as a weakly interacting gas of bosons in two dimensions. Since the 
thermal deBroglie wavelength in two dimensions varies inversely with mass, the extremely 
light mass of these bosonic particles means that the critical temperature for superfluidity 
can in principle be 100 K or above for experimentally achievable number densities. 

In a translationally invariant two-dimensional system, without a trap, superfluidity oc- 
curs via a Kosterlitz-Thouless superfluid (KTS) transition. Experiments on untrapped 
systems^^ have shown promising indications of the onset of spontaneous coherence effects. 
This can be viewed as a type of Bose- Einstein condensation (BEC), with coherence length on 
the order of the size of the cloud of particles, what is sometimes called a "quasicondensate."- 
It is possible, however, to have a true Bose-Einstein condensation (BEC) quantum phase 
transition in two dimensions, if there is a confining potential.- 1 ^ Recently, an experimental 
method has been demonstrated for creating such a confining potential trap in a 2D exciton- 
polariton system, in which the exciton energy is shifted using a stress-induced band-gap 
shift—, and evidence for Bose-Einstein condensation of polaritons has been observed in this 
system^. In these experiments, the trap is macroscopic, about 30 microns across compared 
to a typical interparticle distance of 0.3 microns, and the spring constant is low enough that 
the spacing between the quantized states hu>o in the harmonic potential is small compared 
to ksT, so that the states may be treated as a continuum. The diffusion length of the 
polaritons is comparable to the trap size, so that we may consider them to be in equilibrium 
spatially. 

The properties of polaritons have been studied in several theoretical works. The theor y of 
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ariton dynamics due to polariton-polariton interaction has been developed in Refs 
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16 and 



15l . The crossover between lasing and polariton coherence has been studied in Refs. 
171 . Polariton superfluidity has been predicted 18 as well as spontaneous linear polarization of 
the light emission.— In these previous studies, the coherent polaritonic phases were analyzed 
in the 2D infinite system. 

In this paper we present the theory of the trapped polariton condensate. The paper 
is organized in the following way. In Sec. [Til the effective Hamiltonian of of microcavity 
polaritons in trapping potential is derived. In Sec. II III the number of polaritons in Bose- 
Einstein condensate (BEC) as a function of temperatures is calculated. The superfluid 
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fraction as a function of temperature is also obtained. Finally, in Sec. IIVI we present our 
conclusions. 

II. THE EFFECTIVE HAMILTON! AN OF OF MICROCAVITY POLARITONS 

IN TRAPPING POTENTIAL 

The polaritons are linear superpositions of excitons and photons. The Hamiltonian of 
polaritons is given by Htot = H exc + H p h + H exc - P h, where H exc is an excitonic Hamiltonian, 
H p h is a photonic Hamiltonian, H exc _ p h is a Hamiltonian of exciton-photon interaction. 
Analogous to the case of Bose atoms in a trap,— ^ in the case of a slowly varying external 
potential, we can make the quasiclassical approximation, assuming that the effective exciton 
mass is not a function of r. This quasiclassical approach is valid only if the characteristic 
Pr ^> h, where P is the momentum and r is the radial coordinate in the trap. This is the 
case in the recent experiments-^ 

The Hamiltonian of 2D excitons in the infinite homogeneous system is given by 

tc=E^)fe + i E uj^+jh^tptp,, (i) 

P Z/1 P,P',q 

where b P and bp are excitonic creation and annihilation operators obeying to Bose commu- 
tation relations. In the first term, e ex (P) = E band — E binding + e (P) is the energy dispersion 
of a single exciton in a quantum well, where E^and is the band gap energy energy, Eluding = 
Ryl = /ie-he 4 / (h 2 e) is the binding energy of a 2D exciton (fx e -h = m e mh/(m e + rrih) is 
the reduced excitonic mass, e is the dielectric constant, and e is the charge of an electron), 
and £o(p) = P 2 /(2M), where M = m e + rrth is the mass of an exciton. In the second, 
interaction term, A is the macroscopic quantization area and [7 q is the Fourier transform 
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in the 



of the exciton-exciton pair repulsion potential. As discussed in Refs. 
low-density limit, the excitons can be treated as pure bosons, with an interaction poten- 
tial that includes effects of the underlying fermion nature of the electrons and holes. For 
small wave vectors (q <^ cl^di where cl 2 d — ^ 2 e/(2/i e -he 2 ) is the effective 2D Bohr radius 
of excitons) the pair exciton-exciton repulsion can be approximated as a contact potential 
Uq ~ Uq = U = 6e 2 a 2 D/e- This approximation for the exciton-exciton repulsion is appli- 
cable, because resonantly excited excitons have very small wave vectors.— Another reason 
for this approximation is that the exciton gas is assumed to be very dilute and the average 
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FIG. 1: a) Dashed line: the energy of the bare exciton vs. r, for the parameters of the 
GaAs/AlGaAs structure described in the text. Dotted line: the energy of the cavity photon mode. 
Solid lines: the upper and lower polariton energy which arises from the mixing of the photon and 
exciton modes, b) The photon fraction as a function of r, for two cases: solid line: exciton energy 
resonant with the cavity photon mode at r = 0; dashed line: exciton energy detuned 10 meV below 
the cavity photon. 

distance between excitons r s ~ (7m) _1 / 2 ^> a 2 D, which implies the characteristic momentum 

q ~ r^ 1 <C d2D- A much smaller contribution to the exciton-exciton interaction is also given 

by band-filling saturation effects,— which are neglected here. 

The spatial dependence of the external field V(r) comes about due to the shifting of the 

exciton energy with inhomogeneous stress;— the photon states in the cavity are assumed 

to be unaffected by stress. In this case the band energy E band is replaced by E band {r) = 

-E'feand(O) + V(r). Near the minimum of the exciton energy, V(r) can be approximated as 
1 2 

The Hamiltonian of non-interacting photons in a semiconductor microcavity is given by^: 

H ph = £ph(P)al>a P , (2) 
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where a P and dp are photonic creation and annihilation Bose operators, and £ p h(P) = 
(c/n)\J P 2 + fi 2 ir 2 L c 2 is the cavity photon spectrum (c is the speed of light in vacuum, Lc 
is the length of the cavity, and n = y/e is the effective refractive index). 
The Hamiltonian of harmonic exciton-photon coupling has the form:— 

Hexc- P h = hQ R ^ apfrp + h.c, (3) 
p 

where the exciton-photon coupling energy represented by the Rabi constant TiVtp depends 
on the overlap between the exciton and photon wavefunction and the semiconductor dipole 
moment.— We neglect anharmonic terms for the exciton-photon coupling. 

The linear part of the total Hamiltonian H tot (without the second term on the right-hand 
side of Eq. (CQ)) can be diagonalized by applying unitary transformations and has the form:— 

#o = £ e LP {P)py F + Y, e UP (P)uUp, (4) 
p p 

where p\> and u P are the Bose creation and operators for the lower and upper polaritons, 
respectively; the energy spectra of the low/upper polaritons are 

SLP/UP(P) = £ jdQ±£dS. T \^ £ph{P) _ £ex{P)) 2 + A \ mR \2 } (5) 



which implies a splitting between the upper and lower states at P = of 2Qp, known as 



the Rabi splitting. For the GaAs cavities used in Ref. [11|, this splitting was approximately 
14 meV. The upper and lower polariton energies from (J5]) are plotted in Figure 1(a) as a 
function of r at momentum P = for a value of 7 chosen to give a fit to the experimentally 
measured curvature of the lower polariton branch. The fit implies 7 = 960 eV/cm 2 for the 
bare excitons. 

The excitonic and photonic operators are defined as 15 

bp = Xpp-p — Cpu-p, tip = Cppp + Xpu-p, (6) 

where pp and up are lower and upper polariton Bose operators, respectively, and Xp and 
C P are^ 

1 1 

Xp = — -. , Cp = , , (7) 

, / I I ( Kl R \ 2 /-■ , f e LP (P)-e ph (P) \ 2 

where |Xp| 2 and |Cp| 2 = 1 — \Xp\ 2 represent the exciton and cavity photon fractions in the 
lower polariton.— Figure 1(b) shows the photon fraction at zone center, | C 1 2 , as a function 
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of r. Further from the center, the exciton energy e ex becomes detuned from the cavity 
photon energy, leading the lower polariton to become more photon-like. Because cavity 
photon lifetime is so much shorter than the intrinsic exciton lifetime (~ 2 ps compared to 
~ 100 ps), the polariton lifetime is proportional to the photon fraction. This implies that 
polaritons at higher energy in the trap have shorter lifetime; in other words, there is an 
evaporative cooling effect. As shown in Fig. 1(b), this effect can be magnified by tuning 
the exciton level below the photon level at the center of the trap, so that the polaritons are 
more excitonic there. 

Substituting the polaritonic representation of the excitonic and photonic operators (J6j) 
into the total Hamiltonian Htot, the Hamiltonian of lower polaritons is obtained: 15 

Htot = J2 £ LP( P )PPPP + T7 H ^P.P'.qPp+qPp'-qPPPP', (§) 



2A P,P',q 



where 



£^P,P',q — — Ap +q Xp/Xp/_ q Xp. (9) 

For the slowly changing confinement potential V(r) = —(Ei ian d — Eunding) + {c/n)Ti'KL~ c 1 + 
)pjr 2 (r is the distance between the center of mass of the exciton and the center of the trap), 
the exciton spectrum is given in the effective mass approximation as 

ef}{P) = e ex (P) + V(r) = (c/n)hrrL d l + \r 2 + ^. (10) 

1 /2 

This quasiclassical approximation is valid if P » h/l, where I = (/}/(Mw )) is the size 
of the exciton cloud in ideal exciton gas and u>q = ^J r y/M. 

At small momenta a = 1/2(M~ 1 + (c/n)Lc/hn)P 2 1 and weak confinement 
potential (3 = ■yr 2 /\hfl R \ 1, the single-particle lower polariton spectrum obtained by 
substitution of Eq. ffTUl) into Eq. in linear order with respect to the small parameters a 
and P, is 

£ (P) w {c/n)hnL c l - \htl R \ + Jr 2 + ^(M _1 + (c/n)L c /ti7r)P 2 . (11) 

By substituting Eq. (flO]) into Eq. ((7j) we obtain Xp ps 1/v2. The condition of the validity 
of the quasiclassical approach in Eq. (JT|, Pr ^> h, is also applied here. 

If we measure energy relative to the P = lower polariton energy {c/r^finL^ 1 — l^fi^l? 
the resulting effective Hamiltonian for polaritons in the parabolic trap in P space in the 
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effective mass approximation has the form: 



(pi \ 7t(0) 

where the sum over P and P' is carried out only over P >> h/l (as only in this case the 
quasiclassical approximation used in Eq. (TTOj) is valid), and the effective mass of a polariton 
is given by 

= \{M~ l + (c/n)L c /hn); (13) 

the effective external potential V e s(r) = \V(r) (i.e., r y e s = 7/2), and the effective polariton- 
polariton pair repulsion potential is given by the hard-core contact potential U e $(r — r') = 
t/gjj S(r — r') = \ Uci5(r — r'). Using the experimental parameters for GaAs/AlGaAs structure 



used in Ref. 
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(E p h = 1.60735 eV) we obtain M e fj = 7.8 x 10~ 5 m , where mo is the 
vacuum electron mass, in good agreement with the value of 7 x 10~ 5 mo obtained from direct 
measurement of the effective polariton mass using angle-resolved photon detection, reported 
in Ref. Ill 



III. BOSE-EINSTEIN CONDENSATION AND SUPERFLUIDITY OF 

MICROCAVITY POLARITONS 

In the real space the effective Hamiltonian for trapped polaritons will look exactly like 
the Hamiltonian of weakly-interacting dilute 2D Bose gas in a confinement: 

r ~ / h 2 X7 2 \ ~ /7 (0) r ~ 

H cS = j dv^(v) f -— + Kff(r)J V(r) + -f- j dv^ (v)^ (v)^(v)^(v) , (14) 

where ^^(r) and ^(r) are real space Bose field operators of creation and annihilation of 
polaritons, correspondingly. 

Although Bose-Einstein condensation (BEC) cannot happen in a 2D homogeneous ideal 
gas at non-zero temperature, as discussed in Ref. |8|, in a harmonic trap BEC can occur in two 
dimensions below a critical temperature T c ° given by ksT® = hn' 1 6-y c s N/M e f[, where is 
the Boltzmann constant, and N is the total number of polaritons in a trap. This expression 
for the temperature of BEC is valid if we neglect the polariton-polariton repulsion, i.e., if 
we assume U^g = 0. 
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FIG. 2: The condensate profile no(r) in the trap, for 7 = 960 eV/cm 2 and total density at the center 
n(0) = 1.2 x 10 9 cm -2 . Solid line: T = 0, using the Thomas-Fermi/Bogoliubov approximation. 
Dotted line: T = 15 K, using the theory presented in the text. 



Neglecting the anomalous averages (Ap'ipj and (tjj^xjj^J via the Popov approximation^!, 
implying the system to be very dilute, that is na\ D <C 1, where n = N/(tcR 2 ) is the total 
density of polaritons, and R is the 2D radius of the trap, the self-consistent equation for the 



non-condensate density n'{r) at temperatures k B T ^> h^j 7eff/ 'M e s can be written asr^- 
M eS k B T 



n'(r) 



2-Khr 



log 1 — exp 



u. 



(0)|2„2 



cff 



(15) 



where no is the total density of condensate, and fi = 2U^n—U^no is the chemical potential 
of the system in the Popov approximation.— For the experimental parameters of interest, the 
size of a trap is R ~ 30 fxm, the effective 2D Bohr radius of an exciton is ci2D — 130 A, and 
7 ~ 10 3 eV/cm 2 , which implies the above equation is valid for T ^> 1 K, which is true in all 
of these experiments. Figure 2 plots the spatial profile of the condensate for the experimental 
parameters of the trap. Note that only the states P 3> h/l can be treated as quasiclassical. 
Since the the characteristic momenta for Bose condensate of weakly interacting particles is 

,V2 ,. „ „ „ „ „ f n „ (o) \V 2 



n 



we can apply the 



P = (2M eS U^nj satisfy to the condition h/l < P = (2M eS U ( eS 
quasiclassical approximation. 

The total number of polaritons in the condensate is given by Nq = N — N', where N' 



n — n (r) 



2n J n'(r)rdr is the total number of non-condensate particles. Assuming no 
and solving the self-consistent equation (fl5l) with respect to the non-condensate density 
n'(r), we obtain the dependence of the fraction N (T)/N of the total number of condensate 
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particles on the temperature T. This is plotted in Fig. [3] for the experimental conditions. 

In the thermodynamic limit N — > oo, the Thomas- Fermi approximation that the kinetic 
energy of the system can be neglected has been proved to be valid for BEC in a harmonic 
trap. 29 For small quasimomenta P <C \/2M cS U^n and small temperatures, the energy 
spectrum of the quasiparticles e(P) is given by^- e(P,r) « c s (r)P, where c s (r) is the sound 
velocity in the Popov approximation^ (c s (r) = y U^g no (r, T) / M e g) . 




FIG. 3: Condensate fraction Nq/N as a function of temperature, for three trap spring constants: 
solid line: 760 eV/cm 2 ; dotted line: 860 eV/cm 2 ; dashed line: 960 eV/cm 2 . 

Since the spectrum of the quasiparticles is a linear sound spectrum satisfying the Landau 
criterium of superfluidity, 3 - superfluidity of the polaritons can occur in the trap. Therefore 
at small temperatures there are two components in the trapped gas of polaritons: the normal 
component and the superfluid component. We define the total number of particles in the 
superfluid component s N s = N — N n , where N n is a total number of particles in the normal 
component. We define N n analogously to the procedure applied for definition of the density 
of the normal component in the infinite system n n ,~ using the isotropy of the trapped 
polaritonic gas instead of the translational symmetry for an infinite system. We imagine 
that a "gas of quasiparticles" rotates in the liquid in the plane perpendicular to the axis 
of the trap with some small macroscopic angular velocity v. In this case, the distribution 
function of a gas of quasiparticles can be obtained from the distribution function of a gas at 
rest by substituting for the energy spectrum of the quasiparticles s(P) — Lz/, where L = r x P 
is the angular momentum of the particle. Assuming Pr/h 3> 1, we apply the quasiclassical 
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approximation for the angular momentum: L pa Pr and e(L,r) = r 1 c s (r)L. The total 
angular momentum in a trap per unit of volume L tot (r) is given by 

' hL ; Ln B (e(r,L)-Lv), (16) 



J tot\ 



(27rkr) 2 ' 

where we assume that at small temperatures the quasiparticles are non-interacting, and they 
are described by the Bose-Einstein distribution function n B (e) = (exp[e/(/c B T)] — l) -1 . For 
small angular velocities, tib (e — Lz/) can be expanded with respect to Li/. Then we get 

Assuming that only quasiparticles contribute to the total angular momentum, we define the 
density of the normal component n n (r) by L tot (r) = n n (r)L , where L = M eS rv is the 
angular momentum of one quasiparticle. For the total number of particles in the normal 
component we obtain 

N n = 2n [ R n n (r)rdr = f* J^ff' rdr, (18) 
Jo Jo n c^{n {r))M cS 

where ((z) is the Riemann zeta function (C(3) — 1.202), and the density of the condensate 

no{r) = n — n'(r) (the density of non-condensate polaritons n'(r) can be obtained from 

Eq. (1151) ). The dependence of the fraction N S (T)/N of the total number of polaritons in the 

superfluid component N S (T) = N — N n (T) on the temperature T is presented in Fig. HI The 

superfluid fraction depends only weakly on the spring constant 7, and in the limit 7^0 

approaches the the superfluid density for a 2D translationally invariant system.— 



IV. DISCUSSION 



In conclusion, at low temperature, the Hamiltonian of 2D exciton polaritons in a slowly 
varying external parabolic potential corresponds directly to the case of a weakly interacting 
Bose gas with an effective mass and effective pair interaction in a harmonic potential trap. 
The condensate fraction and the superfluid component are decreasing functions of tempera- 
ture, as expected, and increasing functions of the curvature of the parabolic potential. The 
mixing with the photon states leads to a decreased lifetime for high-energy states, that is, 
an evaporative cooling effect, but does not fundamentally prevent condensation. 

The results given here are comparable to those of the experiments, but do not correspond 
exactly. The condensate peak seen in the experiments with a trap has approximately 15 
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FIG. 4: Superfluid fraction N s /N as a function of temperature for the same three trap spring 
constants as Fig. 3. Dashed-dotted line: the superfluid fraction in the limit 7 — ► 0, namely, the 
translationally invariant 2D case. 

/im full width at half maximum, while the peak shown in Fig. 2 comparable conditions has 
width of approximately 30 /im. The most likely reason is that the mean-field shift due to the 
repulsive interaction between particles U e sn is known to strongly overestimate the actual 
energy shift by as much a factor of ten,— because anticorrelation of the excitons tends to 
reduce the average interaction potential. 

The condensate fraction obtained is a decreasing function of the characteristic potential 



of the interparticle repulsion, which corresponds to the results obtained in Ref. 
authors of Ref. 



. The 



33j showed that at finite number of bosons N, the interparticle repulsion 



suppresses the temperature of BEC, and in the thermodynamic limit N — > 00 the interpar- 
ticle interaction eliminates BEC at finite temperatures. Since we consider the very dilute 
gas of a finite number of polaritons with weak repulsion (weakly-nonideal Bose-gas), the 
increase of the interparticle repulsion results in the increase of the non-condensate fraction 
(Eq. (fT5|) ) at the fixed finite temperature, which agrees with the results of Ref. 33]. 

In our calculations, we have assumed thermal equilibrium. Since the polariton lifetime is 
short, one may question this assumption. The condition for thermal equilibrium, however, 
is simply that the time scale for thermalizing collisions be short compared with the particle 
lifetime. Porras et al.— have shown than the time scale for polariton-exciton scattering 
can be fast enough for a thermalized distribution of polaritons to exist in lowest fc-states. 
Although polaritons have very short lifetime, thermodynamic equilibrium can be achieved 
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in the regime of the strong pump. Polariton- polariton interactions can help overcome of the 
bottleneck and lead to large occupation numbers of the ground state. However, we cannot 
rule out that consideration of pump and decay in a steady state may lead to differences 
from the results presented here. To give an example, the renormalized dispersion of BEC 
of particles with infinite lifetime is Bogoliubov-like, while in the steady state of the system 
with pump and decay it is very different.— This consideration of the influence of the decay 
on the BEC may be a subject of further studies of a trapped gas. 

The spin polarization is important not only for the excitations, but for the condensate 
itself. In Ref. js] the dynamics of the spin of the polariton BEC was analyzed in details. 
It was shown that polariton-polariton interactions lead to the polarization dephasing in 
spatially confined systems. The influence of spin on the phase transitions of the trapped 
polaritonic gases are the subject of further research. 
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